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Abstract. The lPetrova I ll2000l) model to calculate pulse profiles is extended to a variable emission height model to 
make it physically self-consistent. In this context variable means that the emission height is no longer considered 
to be the same for different magnetic field lines. The pulse profiles calculated using this new model seem to be 
less realistic due to a focusing effect and cannot be used to fit (typical) multifrequency pulsar observations. Apart 
from the focusing effect the general morphology of pulse profiles is not greatly affected by introducing a variable 
emission height. Additional extensions of the model will be needed to be able to fit observations, and several 
suggestions are made. 
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1. Introduction 



lArons fc Barnard! l)l986|) have derived the dispersion re- 
lation for three wave modes which can propagate through 
the plasma of a pulsar magnetosphere: the ordinary sublu- 
minous mode (subluminous 0-mode), the ordinary super- 
luminous mode (superluminous 0-mode) and the extraor- 
dinary mode (X-mode). The X-mode does not suffer re- 
fraction, but refraction of the sublumin ous 0-mode can be 
considerable in pulsar magnetospheres l|Barnard fc Arons I 
Il986l|) . The subluminous 0-mode cannot escape the pul- 



sar magnetosphere due to Landau damping, s o it does not 
contribute directly to the observed emission. iLvubarskii I 
l)l996|) has shown that the subluminous 0-mode can be 
converted into the superluminous 0-mode - which can es- 
cape the magnetosphere - by induced scattering; off pla sma 
particles. As pointed out bv lBarnard fc Arons I l)l986|l re- 
fraction of the superluminous 0-mode is less severe than 
for the subluminous 0-mode. It can, however, be impor- 
tant in the presence of a transverse plasma density gradi- 
ent. 

For the superluminous 0-mode IPetrova' ] ll2n0(t) (here- 
after P2000) shows how pulse profiles can be calculated 
taking into account the transverse plasma density gradi- 
ent. This model demonstrated that complex profiles can 
be produced by a "simple" ring shaped emission region 
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(as predicted bv lRuderman fc Sutherland Il975|) . and thus 
that the wealth of observed pulse profile shapes may be 
due to different magnetospheric conditions rather than 
more complex emission region shapes. Furthermore it was 
shown that the observed phenomenon of high frequency 
core splitting could be an effect of refraction. 

The emission height is an important ingredient in 
calculating pulse profiles. The emission height is fre- 
quen cy depeiident" i.e. there is radius-to-frequency map- 
ping l|Cordes 119781) . Plasma waves with higher frequencies 
are excited closer to the star. The observed frequency de- 
pendence of pulse profiles is often very complex, perhaps 
more complex than can be expected from just radius-to- 
frequency mapping. Because refraction itself is a frequency 
dependent phenomenon, a more complex frequency de- 
pendence of pulse profiles can be expected if refraction 
is important in pulsar magnetospheres. Other effects that 
can be understood by taking into account refra ction are 
the o ccurrence of orthogonal polarization modes l|Petrova I 
l200l|) and the spectral breaks of pulsars ijPetrova Il2002(l . 

To link the observed pulse profiles to the shape of the 
emission region, so as to be able to check emission theories, 
one must know the refractive properties of pulsar mag- 
netospheres. This calls for the development of improved 
refraction models. As noted by P2000, the emission sur- 
face at one observing frequency should be, strictly speak- 
ing, an isodensity surface of the plasma distribution. Yet, 
for simplicity, a constant emission height (CEH) was as- 
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Magnetic Axis 




Fig. 1. The ray at position (r, x) is propagating in the 
direction 9. This ray was emitted at {ro,Xo) and the field 
hne through this point is indicated by Xs- The plasma 
density peaks at the characteristic field lines indicated by 
Xc- The angles Xc and Xs are defined at r = 1 and the 
angles Xo and at the emission height tq. 

sumed in P2000, in the expectation that the qualitative 
features of profile formation would not be sensitive to that 
assumption. In the present paper we do adopt a surface of 
constant density as required for self-consistency of the re- 
fractive model, and we investigate the effects on the pulse 
morphology. This 'variable' emission height (VEH) ap- 
pears to introduce a focusing effect which causes the pro- 
files to have unrealistically sharp edges. As a consequence, 
the VEH model cannot be used to fit multifrequency pul- 
sar observations without relaxing additional restrictive as- 
sumptions, a number of which are discussed at the end of 
the paper. 



2. Refraction model 

2.1. The ray equations 

The refraction model below is essentially that of P2000, 
and we refer to that paper for details. The plasma distri- 
bution and the magnetic field are assumed to be axisym- 
metric around the magnetic pole, so the refraction model 
can be described in two dimensions. A position on a ray 
trajectory is indicated by the polar coordinates r and x, 
and the direction along the trajectory by 9 (see Fig. 

The geometrical optics description applies and the 
time evolution of these quantities is given by the 



Hamilton equations. For a highly magnetized ultrarel- 
ativistic electron-positron plasma, which is cold in the 
pr oper restframe, the disper sion relation has been derived 
bv lArons fc Barnard Nl98(il) and the a ssociated Hamilton 
equations bv iBarnard fc Arons I l|l98(jl . On the condition 
that the plasma flows with the same velocity for all field 
lines and when rays are emitted parallel to the local (dipo- 
lar) magnetic field, the dispersion law describing the two 
ordinary wave modes can be written as (P2000) 
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The radial plasma density derivative has been omit- 
ted, because it can be neglected for the plasma density 
we will adopt (P2000). The parameter 77 is related to the 
component of the wave vector k in the direction of the 
local magnetic field and is defined as 



V = 27^(1 - "11) 



(3) 



with 7^1 the Lorentz factor of the outflowing plasma, 
and n|| — ck^^/uj where uj is the frequency of the plasma 
wave. The refractive index is n = (riy +n'^y^'^ , where par- 
allel/perpendicular is with respect to the local magnetic 
field. It is assumed that ny is such that 77 ^ 27^. The 
plasma waves are assumed to be generated close to the 
local Lorentz-shifted plasma frequency u)py/j, 



m 



(4) 



with e the electron charge, m the electron rest mass and 
Np the particle number density (electrons plus positrons) 
of the plasma. The ratio 



/ 



(5) 



should then be close to unity. 

Eqs. and (|2Jl are normalized, i.e. the coordinates r, 
Xn and 0,1, as well as the plasma number density distri- 
bution N, arc normalized to their values at the emission 
height (so x — XoXn and 9 — 9o9n)- The emission height 
can be different for different rays as will be discussed in 
Sect. 2.3, so in this context the emission height is the emis- 
sion height of the particular ray that is being considered. 
The values of / and x at the emission height are denoted 
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Fig. 2. Rays (dashed) with the observing frequency are generated at the emission height and are refracted toward 
lower plasma densities. The scale is in units emission height and the gray scale indicates the hollow cone plasma 
density. The plasma density is proportional to the numbers of the gray scale bar. The plots on the left are for the 
CEH model, the plots on the right are for the VEH model, the top plots are calculated with Xc = 0.03 and the bottom 
plots with Xc = 0.01. The other input parameters are: 7 = 30, /o = 0.5, ei = 3, £2 = 4. 



as /o and xo respectively. All angles are assumed to be 
small compared to 1 throughout this paper, so the prop- 
agation direction of the plasma waves should always be 
nearly parallel to the magnetic axis. 

For the superluminous 0-mode which is considered 
here, n|| cannot be larger than 1, therefore rj is required 
to be positive (or zero). At the emission height (where 
On — Xn — N — 1) the solution of the dispersion relation 
follows immediately and we have for the superluminous 
0-mode 



_ r 2//0 - 1 for < /o < 2 
10 for /o>2 



(6) 



Note that 770 is continuous at /o = 2. The solution of the 
dispersion relation applicable above the emission height is 
given by the general solution of the cubic ljA.l|l . 

Eqs. Q and ^ describe, to first order in Xn and 6'„, 
the refraction of an ordinary (both the sub- and super- 
luminous) plasma wave. The two differential equations 
for Xnir) and Onir) can be solved numerically if 77 (r) 
is known. As noted above, 77 can be calculated analyti- 
cally from the dispersion equation. We use a fourth or- 
der Run gc-Kutta method with adaptive stepsizc control 
l)Press et al.. .198fii) to solve the set of differential equa- 




tions. For the plasma density distribution we adopt the 
hollow cone model (P2000) 



(7) 



where Ni, is the particle number density at (r = 1, x = 
Xc)- This plasma density is Gaussian shaped around a 
"characteristic field line" indicated by Xc- The decrease 
of the plasma density may be different for the inner and 
outer regions, so we set 



ei for Ixl < XcVr 
£2 for Ixl > XcVr 



(8) 



Eq. ^ do not contain ro , so the whole problem is indepen- 
dent of the scaling of r. However, in Eq. lO Xc is defined 
at r = 1. The normalized plasma density is given by 
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and its derivative with respect to Xn is 
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Fig. 3. The final propagation direction Of versus Xs for botli tlie VEH and tlie CEH model, for two values of Xc- In 
these plots Xs ranges from Xc/5 to 5xc/2 and the other parameters are the same as for Fig. |21 The final propagation 
direction 0f is in radians and Xs in units Xc- 



The parameters required to calculate a single ray trajec- 
tory are those of the plasma density distribution (xc £i 
and £2), the plasma outflow Lorentz factor 7, the fre- 
quency of the plasma wave (expressed in /o) and the start 
position xo of the ray. Solving Eqs. ^ and (O with the 
start condition Xn = = 1 will give the ray trajectory. 

Plasma waves are refracted toward lower plasma den- 
sities in the magnetosphere until refraction becomes in- 
efficient due to the decreasing plasma density along its 
trajectory. As can be seen in Fig. |21 refraction results in 
a redistribution of rays; i.e. the rays are no longer equi- 
spaced above the emission height and two "conal compo- 
nents" of outer rays and a "core component" of inner rays 
are formed. 

2.2. Calculation of the pulse profiles 

The effect of refraction is quantified in a plot of the final 
ray direction (Fig.O. Here the propagation direction 6*/ at 
a height where refraction has become inefficient is plotted 
versus Xs- The initial ray position (ro, xo) corresponds to a 
value of Xs by tracing back the field line from the emission 
height to r = 1 (see Fig.^. 

Both inside and outside the characteristic field line 
cone refraction is toward lower plasma densities, which 
results in a steepening of the final ray direction plot near 
Xs = Xc- Inside the plasma cone the rays are refracted to- 
ward the magnetic axis and the innermost rays may even 
intersect the magnetic axis; in that case the final ray di- 
rection plot crosses the line Of = 0. 

For small values of 0/, rays originating from several 
discrete values of Xs leave the magnetosphere in the same 
direction and at the corresponding pulse longitude differ- 
ent parts of the emission ring can be observed simulta- 
neously. Note that the final ray direction curve for the 
opposite half of the emission ring is found by mirroring 
the curve with respect to the line Of =0. If the final ray 
direction plot crosses this line, some parts of both sides 
of the emission ring have the same df.ln that case both 
sides of the emission ring can be observed simultaneously, 
if the impact angle (3 is small enough. 

If the curve in Fig.|21is horizontal at the Of value cor- 
responding to the fine of sight (6'los), a large part of the 



emission surface is observed simultaneously while if the 
curve is steep only a small part is observed. This means 
that the observed intensity in the pulse profile is propor- 
tional to the value of dxs/dOf at Of = 0los which is just 
an energy conservation argument (P2000). 

Apart from refraction effects the pulse profile will de- 
pend on the intensity distribution at the emission height. 
If the pair production is someh ow related to the observed 
coherent microwave radiation ijRuderman fc Sutherland I 
1 975i) , then similar distributions for the plasma density 
and the intensity at the emission height can be expected 
such as (P2000) 

W'.o = exp (^-£T (^^i^^^yy (11) 

This corresponds to an emission ring which peaks at the 
characteristic field lines and its thickness is set by T. For 
T = 1 the intensity distribution follows exactly the plasma 
density distribution. The shape is Gaussian as a function 
of the field line parameter Xs, but the choice of another 
parameter (such as the length along the emission surface) 
is also conceivable. However for simplicity the parameter 
Xs is used. As will be discussed later on, the conclusions 
do not depend on this choice. 

The refraction model is axisymmetric around the mag- 
netic axis, so the beam-pattern of the pulsar is also ax- 
isymmetric around the same axis. The shape of observed 
pulse profiles depends on how the line of sight cuts the 
magnetic pole of the star. We only consider the most sim- 
ple geometry; i.e. the magnetic axis is orthogonal to the 
rotation axis (a = 90°) and the line of sight cuts the 
magnetic pole centrally (impact angle /3 — 90°). For this 
geometry the pulse longitude is equal to the final ray 
direction Of. Because of this choice of geometry and the 
axisymmetry, all the information of the beam-pattern is 
in the calculated pulse profiles. The model itself is inde- 
pendent of a and /3, only the mapping between (j) and Of 
changes. 

2.3. Variable emission height model 

The model described above may be applied with both a 
constant (CEH) and a variable emission height (VEH), 



p. Weltevrede et al.: The effect of a variable emission height on pulse morphology 
Constant emission height (CEH) Varable emission height (VEH) 

X =0.02 z,=ooi Zj.=0.0? Z,. =0.02 x^=o.oi 






10 -10-K -6 -4 -2 2 4 [> S 10 -10 -S -6^-2 2 




Fig. 4. Pulse profiles for different Xc (a small Xc corresponds to a high observing frequency) for both the CEH and 
the VEH model. For the top row Xs ranges from Xc/2 to 3xc/2 with T = (all field lines having equal intensity) and 
for the bottom row Xs ranges from Xc/5 to 5xc/2 with T = 1 (intensity coupled to the plasma density). The other 
parameters are the same as in Figs. El and 13 The scale is such that the integrated intensity is the same for all profiles. 



but (as we will argue) a VEH is needed to make the model 
self-consistent. The requirement of a VEH was not met in 
P2000; it is the basic conceptual difference between the 
model presented here and the P2000 model. Its effect on 
the pulse profiles turns out to be appreciable, as discussed 
below. 

The emission height can be derived when a plasma den- 
sity distribution has been specified. The plasma density 
decreases as ^ so the local plasma frequency decreases 
away from the star resulting in the excitation of plasma 
waves with higher frequencies closer to the star. This re- 
sults in rays propagating in a direction which is more 
aligned with the magnetic axis at the emission height. But 
there is another effect involved in the frequency depen- 
dence of the pulse profile morphology: refraction becomes 
more prominent. 

The assumption that both /o and 7 are constant im- 
plies that plasma waves of one particular frequency are 
generated at one particular equi-plasma density surface 
(Eqs.0]and|Sll. Because a transverse plasma density gradi- 
ent is needed for refraction, the emission height of a given 
frequency varies with polar angle xo ■ If the magnetic field 
is dipolar, we have 

Xo = Xs^/ro, (12) 

where Xs is shown in Fig.^ Combining the plasma density 
(Eq. I?!) with Eqs. 10} and lO at the emission height, and 
using Eq. ifT^ leads to the following expression for the 
emission height 



R exp — £ 



Because 7 and /o are assumed to be constant, R should 
be constant and vq is not constant. As noted earlier, the 
whole problem is independent of the scale of r, so R can 
be set equal to 1. The frequency dependence of the pulse 
profiles is then in the parameter Xc (P2000) 



Xc oc 



-1/3 



(14) 



(13) 



The ray trajectory is solved as a function of the distance 
to the star, expressed in units of the emission height and 
different Xc correspond to relative observing frequencies. 
The physical emission height can be calculated from Eq. 
(|13|l when Ni, and uj are specified. 

The emission surface specified by Eq. (|13() corresponds 
to an isodensity surface, so the plasma density distribution 
has a more prominent role in this VEH model than in 
the CEH model. Apart from causing refraction, it also 
determines the shape of the emission surface. 

Refraction becomes more severe for the inner and outer 
rays in the VEH model, because the plasma gradients are 
larger at lower emission heights. Moreover a lower emis- 
sion height implies that the rays are emitted closer to, 
and are initially propagating more aligned with the mag- 
netic axis. For the inner rays this means that the rays 
can intersect the magnetic axis more easily. For the outer 
rays there are two counteracting effects. A lower emission 
height implies that the rays are refracted in a more out- 
ward direction, but at the same time the rays are also 
emitted more aligned to the magnetic axis. 



3. Results 

Model calculations of pulse profiles for both a VEH and 
a CEH are presented in Fig. ^ for the most simple geom- 
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etry (/3 = 0° and a — 90°). For this geometry the pulse 
longitude cj) is equal to the final ray direction 9f. 

Observationally the core component behaves differ- 
ently from conal components, both in the frequency de- 
pendence of its morphology and in its polarization proper- 
ties. This is what can be expected from refraction (P2000), 
because the core component consists of "mixed" rays; i.e. 
the order of the beams changes. This is true for both a 
VEH and a CEH. 

In Fig. 121 one can see that refraction becomes more 
prominent for higher frequencies (lower Xc)- This can also 
be seen in Fig. |31 where the final propagation direction of 
rays versus Xs is plotted. The curve becomes more complex 
for lower Xc- Besides this refractive effect, a lower emission 
height implies smaller propagation angles 9 at the emission 
height, resulting in narrower pulse profiles with increasing 
frequency (decreasing Xc) in Fig. ^ This is again true for 
both a VEH and a CEH. 

The pulse profiles in Fig. 01 for a CEH are more spiky 
than the pulse profiles presented in P2000. The main rea- 
son for this is the higher resolution of the calculations 
presented here. 

There are three reasons why the profiles, for both a 
VEH and a CEH, are spiky. First of all, if the intensity 
distribution at the emission height is flat, the emission 
ring has sharp edges resulting in sharp edges in the pulse 
profile. This effect can be reduced by making T larger (see 
Eq. as is seen in the bottom row of Fig. 0] Making 
the parameter T larger results in the edge of the intensity 
distribution becoming Gaussian blurred. 

The second effect is caused by rays crossing the mag- 
netic axis, so this applies especially to the core component. 
This means that at certain pulse longitude both sides of 
the emission ring can be seen simultaneously. When the 
number of sides visible changes at a particular pulse longi- 
tude, a step in intensity appears in the pulse profile. This 
effect can again be reduced by increasing T as can be seen 
in the bottom row of ^ 

The last effect contributing to the spikiness of the pro- 
files is a focusing effect. If a large patch of the emission 
ring is focused at one pulse longitude, a peak is observed. 
This focusing effect corresponds to a horizontal part in 
the final ray direction curve; rays emitted at a range of Xs 
are focused to a single Of.At high frequencies this can be 
seen for the innermost rays (Fig.j^l, causing peaks in the 
core component. 

By comparing the pulse profiles for the VEH and the 
CEH model in Fig. ^ the most striking difference is the 
conal components. The edges of the VEH profiles are very 
sharp, and introducing a large T will not reduce their 
sharpness. The reason can be found in Fig. [31 The curves 
for the VEH model show a global maximum at Xs ~ 1-5. 
Because it is a maximum, there is focusing and because 
the maximum is global, the peaks occur at the edges of 
the profile. 

As discussed in Sect. 12.31 there are two counteract- 
ing effects for the outer rays in the VEH model. A lower 
emission height makes the outward directed refraction 
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Fig. 5. Pulse profiles calculated for the case of no refrac- 
tion. For the top row Xs ranges from Xc/2 to 3xc/2 with 
T = and for the bottom row Xs ranges from Xc/5 to 
5xc/2 with T = 1. These profiles have been calculated for 
Xc = 0.02. The other the parameters and the normaliza- 
tion of the profiles are the same as in Fig. ^ 



stronger, but the propagation angle is more aligned with 
the magnetic axis at the emission height. The reason for 
the global maximum is that the latter effect dominates for 
the outermost rays. 

This focusing effect due to the variable emission height 
is also visible in the pulse profiles of Fig. which were cal- 
culated without using refraction. The focusing is caused 
by the geometry of the emission surface, not by the in- 
tensity distribution at the emission height (Eq. . This 
means that this focusing is independent of the precise form 
of Eq. Hll|) , and therefore also of the choice to use the field 
line parameter Xs instead of for example the length along 
the emission surface in this equation. 

The edges of the profiles are produced by rays emitted 
from Xs ~ 1-5 and the rays are focused, so there should 
be only very little radiation produced at Xs ~ 1.5 to avoid 
the sharp edge. This means that the emission ring should 
be very thin, so T should be large. At high frequencies 
T should be at least « 5 and at the lowest frequencies 
(Xc = 0.03) T should be at least sa 15. A large T phys- 
ically means that only the middle part of the emission 
ring is producing coherent microwave radiation, although 
the whole ring is producing streams of particles. Such a 
scenario is in conflict with the expectation that pair pro- 
duction and coherent emission are related. 
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A VEH leads to stronger refraction. Besides introduc- 
ing a VEH refraction can also be increased by changing the 
values of El, e^^ fo or 7 in the CEH model. Experimenting 
with a range of values of these parameters did not lead 
to the formation of the sharp edge of the profiles with a 
CEH. Therefore the focusing effect is a typical property 
of the VEH model. 



4. Discussion 

Contrary to the expectation expressed in P2000, the qual- 
itative features of profile formation turn out to be different 
for the VEH and the CEH refraction models. Although the 
VEH is a physical improvement in the sense that it makes 
the emission model self-consistent, the profiles obtained 
are less realistic. The model, therefore, needs further im- 
provements before it can serve as a tool to fit (typical) 
multifrequency pulsar observations. 

The most pronounced difference between the CEH and 
the VEH model is that for the VEH model the rays emit- 
ted at the outside of the emission surface do not form the 
edges of the pulse profile. The edges of the pulse profiles 
in the VEH model are generated by a focusing effect caus- 
ing the edges to be sharp. If the thickness of the emission 
ring at the emission height were much thinner than the 
thickness of the plasma cone at the emission height the 
sharp edges would disappear, but this seems physically 
unrealistic. 

It must be noted that the results depend strongly on 
the plasma distribution adopted. The density profile not 
only causes refraction, but it also determines the shape of 
the emission surface. If the plasma density falls off more 
slowly than the Gaussian distribution assumed here, the 
results may be more realistic although in that case refrac- 
tion will be less prominent. 

Several other effects could contribute to smoother 
pulse profiles. There are probably more frequencies gen- 
erated at one point in the magnetosphere, so there would 
be a fo range rather than a fixed fo value. Also the rays 
are not emitted strictly aligned with the magnetic field 
lines, rather there will be an elementary beam pattern 
of finite angular width. A beam pattern with a width of 
can be considerable compared with the pulse width 
(for a plasma outflow Lorentz factor 7 w 30 the beam is 
about 2° wide). If the outflow Lorentz factor is different 
for different field lines, the shape of the emission surface 
is changed. Moreover refraction becomes more complex, 
because it depends on gradients of the 7 factor as well a s 
gradients in the plasma density ijBarnard fc Arons1ll986() . 
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cubic. The solution for the superluminous 0-mode (the 
solution with positive 77 at the emission height) is 

77 = s+ + s_-^, (A.l) 

where 

/ \ 1/3 

s± = ^r±^/q^ + r'^j 

n - £i _ iil (A.2) 

y - 3 9 3 

r = |(aia2 - 3ao) - |f 
and 

a2 = 2 + flo 
N 

ai = l-4--2+2ao (A.3) 
Jo 

9 2 2//1 \2 
flO = "4X07 (6'n~Xn) • 

Acknowledgements. We thank Svetlana Petrova for making 
available her code and for her valuable comments as referee 
of this article. We also thank her as well as Joeri van Leeuwen, 
Ramachandran and John Barnard for constructive discussions. 

References 

Arons, J. & Barnard, J. J. 1986, ApJ, 302, 120 
Barnard, J. J. & Arons, J. 1986, ApJ, 302, 138 
Cordes, J. M. 1978, ApJ, 222, 1006 
Lyubarskii, Y. E. 1996, A&A, 308, 809 
Petrova, S. A. 2000, A&A, 360, 592 
Petrova, S. A. 2001, A&A, 378, 883 
Petrova, S. A. 2002, A&A, 383, 1067 

Press, W. H., Flannery, B. P., Teukolsky, S. A., & Vetterling, 
W. T. 1986, Numerical Recipes: The Art of Scientific 
Computing (Cambridge: Cambridge University Press) 

Ruderman, M. A. & Sutherland, P. G. 1975, ApJ, 196, 51 



Appendix A: Analytical solution of the normalized 
dispersion relation 

The normalized dispersion relation is of the third de- 
gree in 77, so the analytical solution of 77 is given by the 



